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Based on recent first-principles computations in perovskite compounds, especially BaTiOs, we 
examine the significance of the Born efi'ective charge concept and contrast it with other atomic 
charge definitions, either static (MuUiken, Bader...) or dynamical (Callen, Szigeti...). It is shown 
that static and dynamical charges are not driven by the same underlying parameters. A unified 
treatment of dynamical charges in periodic solids and large clusters is proposed. The origin of the 
difference between static and dynamical charges is discussed in terms of local polarizability and 
I delocalized transfers of charge : local models succeed in reproducing anomalous effective charges 

thanks to large atomic polarizabilities but, in ABO3 compounds, ab initio calculations favor the 
physical picture based upon transfer of charges. Various results concerning barium and strontium 
titanates are presented. The origin of anomalous Born effective charges is discussed thanks to a 
band-by-band decomposition which allows to identify the displacement of the Wannier center of 
separated bands induced by an atomic displacement. The sensitivity of the Born effective charges to 
microscopic and macroscopic strains is examined. Finally, we estimate the spontaneous polarization 
in the four phases of barium titanate. 

> 

■ I. INTRODUCTION 

I For a long time, there has been a continuing interest in the definition of atomic charges in solid state physics as 

■ well as in chemistry [Q-Q . This interest lies essentially in the fact that the concept of atomic charge naturally arises 
00 \ in a large diversity of frameworks and is frequently helpful for a simple description of solids and molecules. The 

variety of contexts in which the charge is involved (IR spectrum analysis, XPS chemical shifts analysis, theory of ionic 
conductivity of oxides, determination of electrostatic potential, definition of oxidation states...) underlines its central 
role but also reveals a concomitant problem: inspired by various models or by the description of various physical 
C I phenomena, many different definitions have been proposed that, unfortunately, are not equivalent 

■ Following a distinction already made by Cochran it seems possible to classify the different concepts into static 
^ \ and dynamical charges. The static charge is an intuitive concept, usually based on a partitioning of the ground-state 
O ■ electronic density into contributions attributed to the different atoms. It is an ill-defined quantity that depends on 

the convention artificially chosen to affect a given electron to a particular ion . On the other hand, the dynamical 
charge is directly related to the change of polarization (or dipole moment, for molecules) created by an atomic 
displacement. This change of polarization is a quantity that can be experimentally measured, at least in principles, 
giving the dynamical charge a well-defined character. 
Ci In order to clarify the concept of atomic charge, it was important to compare on practical examples the numerical 

results obtained from its different definitions. Recent studies of the statistical correlation between various atomic 
charges using a principal component analysis have suggested that the different definitions arc not independent but 
correspond to different scales driven by a unique underlying physical factor jdj. If this assertion seems plausible as 
far as static charges are concerned, we will argue that the dynamical charge should not reduce to the same physical 
factor but should also depend on an additional parameter: the rate of transfer of charge, influenced by the bonding 
with the other atoms of the system and additionally, for large systems, by the condition imposed on the macroscopic 
electric field. 

The Born effective charge tensor Z*^-^) - alias transverse charge -, that is at the center of the present study, is 
a dynamical quantity introduced by Born |^ in 1933. In solid state physics, it is since a long time considered as a 
fundamental quantity because it monitors the long-range Coulomb interaction responsible of the splitting between 
transverse and longitudinal optic phonon modes During the seventies, the Born effective charges were already 
investigated and discussed within empirical approaches (see for example Harrison [^). More recently, they became 
accessible to first-principles calculations [|^-^ , and accurate values have been reported for a large variety of materials. 

For the case of ABO3 compounds, old experimental data and empirical studies had suggested that the 
amplitude of the Born effective charges should deviate substantially from the amplitude of the static atomic charge. 
Surprisingly, this result remained in the dark until first- principles calculations confirmed that the components of Z*^^^ 
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are anomalously large in these oxides |ri|-[i3[|. It was observed that the components of Z*^'^'> can reach twice that 
of the nominal ionic charges. This result reopened the discussion on the physics of the Born effective charges and 
different recent studies tried to clarify the microscopic processes monitoring the amplitude of Z*'^'^\ 

In this paper, we first clarify the relationship between various atomic charges. We then present results concerning 
BaTiOa and SrTiOa in order to illustrate how a careful analysis of the Born effective charges can teach us interesting 
physics concerning these compounds. It reveals the mixed ionic and covalcnt character of the bond |14 1^. It 
allows to visualize the mechanism of polarization as electronic currents produced by dynamical changes of orbital 
hybridizations |^,^. It also clarifies the origin of the giant destabilizing dipole-dipole interaction producing the 
ferroelectric instability of these materials . 

In Section II and III, we contrast the concepts of static and dynamical charges and we reintroduce the Born 
effective charge that is at the center of the present discussion. In Section IV, we compare various results obtained 
within different frameworks for the cubic phase of BaTiOa and SrTiOs. The bond orbital model of Harrison is 
explicitely applied to SrTiOa (Appendix A) . We also discuss the origin of the large anomalous contributions in terms 
of local electronic polarizability and dynamical changes of orbital hybridization. A decomposition of the role played 
by the different bands is reported in Section V. Section VI is devoted to the evolution of the Born effective charges 
in the three ferroelectric phases of BaTiOs as well as in the cubic phase under hydrostatic pressure. This points out 
the role of the anisotropy of the atomic environment on the amplitude of Z*^'^\ Finally, in Section VII, we report 
the evolution of the effective charges all along the path of atomic displacements from the cubic to the rhombohedral 
phase and we estimate the spontaneous polarization of the three ferroelectric phase of BaTiOs. 



II. THE CONCEPT OF STATIC CHARGE 



Intuitively, the atomic charge may first appear as a static concept. The charge associated to an isolated atom is a 
well-defined quantity. The purpose of defining static atomic charges is therefore to extend this notion to molecules 
and solids. For these cases, the challenge basically consists to replace the delocalized electronic density by localized 
point charges associated to each atom. This could a priori be performed from electronic density maps obtained 
experimentally or theoretically. However, as already mentioned by MuUiken Q in 1935, "there are some difficulties 
of giving exact definition without arbitrariness for any atomic property". During the seventies, Cochran similarly 
emphasized that the partition of the electronic distribution into atomic charges can only be done unambiguously 
when "boundary can be drawn between the ions so as to pass through regions in which the electron density is small 
compared with the reciprocal of the volume inclosed" . This is never the case in practice, and especially when there 
is appreciable covalent bonding. For most of the solids and molecules, there is consequently no absolute criterion to 
define the static atomic charge and a large variety of distinct definitions have been proposed that are not necessarily 
quantatively equivalent (see for instance Ref. 1^,^). 

As an illustration, different approaches have been considered in order to evaluate the amplitude of the static atomic 
charges of barium titanate. Some results are summarized in Table ^, where different atomic charges are reported in 
comparison with the nominal charges expected in a purely ionic material (+2 for Ba, +4 for Ti, ~2 for O). Some 
of them were obtained from empirical models; others were deduced from first-principles. The static atomic charges 
of Ref. |6| were deduced by Harrison within his bond orbital model using universal parameters and neglecting the 
interactions with the Ba atom. The atomic charges reported by Hewat jl^ were approximated from a model initially 
designed by Cowley for SrTiOa. The charges reported by Khatib et al. |jl^ have been obtained in a shell- model 
context. In two references, Turik and Khasabov ||2^,|2^ estimated the charges from the Madelung constant thanks 
to a fit of the crystal energy with shell-model parameters. Michel-Calendini et al. proposed charges from a 
population analysis of the Xa electronic distribution of a TiOg cluster, assuming a charge of +2 on Ba. Cohen and 
Krakauer |23) deduced the atomic charges from a fit of the DFT-LDA electronic distribution by that of overlapping 
spherical ions (generated according to the potential induced breathing model) for different ionic configurations. Xu 
et al. reported values deduced from a MuUiken population analysis of a self-consistent OLCAO calculation [ p5| . 
In another reference p6[ | , Xu et al. proposed different values by integrating the electronic charges in spheres centered 
on the ions, and partitioning rather arbitrarily the remaining charge outside the spheres following a method proposed 
in Ref. |2|,|2|. 

The results of Table ^ are not quantitatively identical and illustrate that there is no formal equivalence between the 
different procedures used to define the atomic charge. However, in agreement with an analysis reported by Meister and 
Schwartz Q for the case of molecules, we observe that the values of Table | have some common features, suggesting 
that the different charges are not independent but should correspond to different scales driven by a common factor. 
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In particular, all the calculations reveal that the charge transfer from Ti to O is not complete. If BaTiOa was a 
purely ionic crystal, the 3d and 4s electrons of Ti would be entirely transferred to the oxygen atoms, yielding a charge 
of +4 on titanium. However, due for instance to the partial hybridization between O 2p and Ti 3d states [p9|-^, 
these electrons remain partly delocalized on the Ti atom so that the static charges on the Ti and O atoms are smaller 
than they would be in a purely ionic material. This delocalization is illustrated in Figure |l|, where we have plotted 
the partial electronic density associated to the O 2p bands. For the Ba atom, the situation is not so clear than for 
titanium but most of the studies suggest similarly that the 6s electrons are not fully transferred to the oxygen. 

From the previous examples, it is clear that, strictly speaking, the static charge does not give a quantitative infor- 
mation. In the study of mixed ionic-covalent compounds, it remains however a useful concept to discuss qualitatively 
the transfer of charges from one atom to the other. As a general rule, the partial covalence reduces the amplitude 
of the static charge. Comparison of a specific static charge in the different phases of a given material or in different 
compounds can therefore give a relevant information on the evolution of the chemical bond . 



III. THE CONCEPT OF DYNAMICAL CHARGE 

As emphasized by Harrison ||^ , "whenever an ambiguity arises about the definition of a concept such as the atomic 
charge, it can be removed by discussing only quantities that can be experimentally determined, at least in principles" . 
If there are some ambiguities to determine the charge to be associated to a given atom, the charge carried by this atom 
when it is displaced is directly accessible from the induced change of polarization (or dipole moment for molecules). As 
it is now discussed, the dynamical charges are defined by the change of polarization induced by an atomic displacement; 
from the viewpoint of Harrison, they appear therefore as more fondamental quantities. 



A. Role of the macroscopic electric field 

In molecules, the change of dipole moment in direction /3 (pp) linearly induced by a small displacement of atom k 
in direction a (tk,q.) is uniquely defined. The proportionality coefficient between the dipole moment and the atomic 
displacement has the dimentionality of a charge and is usually referred to as the atomic polar tensor (APT) : 



dpi3 

This concept was introduced by Biarge, Herranz and Morcillo [^6[-p8| for the interpretation of infra-red intensities 
measurements. Later, Cioslowski [ p9|j40| introduced a scalar charge in connection with this tensor : it is the generalized 
atomic polar tensor (GAPT) defined as one-third of the trace of the atomic polar tensor. 

In periodic systems, equivalent atoms appear in the different unit cells and the definition of the charge can be 
generalized. A dynamical charge tensor is conventionally defined as the coefficient of proportionality at the linear 
order between the macroscopic polarization per unit cell created in direction [3 and a rigid displacement of the sublattice 
of atoms K in direction a, times the unit cell volume fio ■ 



Zl.aP - ^0 ^ (2) 



We note that ^o-'P can be interpreted as a dipole moment per unit cell. As one k atom is displaced in each unit cell, 
in the linear regime, this definition is equivalent to Eq. (|l|) : it corresponds to the change of dipole moment induced 
by an isolated atomic displacement. However, contrary to the case of molecules, in macroscopic systems, the previous 
quantity is not uniquely defined. Indeed, the change of polarization is also dependent on the boundary conditions 
fixing the macroscopic electric field £ throughout the sample. Basically, we can write : 
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As the electrostatics imposes a relationship between macroscopic polarization, electric and displacement fields 



V^^£^ + inVa = V (4) 
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we can deduce the following equivalent expression : 



^K,al3 — ^0 



£=0 



in 



(5) 



Depending on the condition imposed on the macroscopic electric field, different concepts have historically been intro- 
duced |6). 

The Born effective charge ^ - alias transverse charge, Z*'^^ - refers to the change of polarization that would be 
observed under the condition of zero macroscopic electric field, so that the second term appearing in the previous 
equation vanishes : 



„*(T) _ ^ dVfj 



(6) 
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The Callen charge - alias longitudinal charge, Z*^^^ - is defined from the change of polarization under the 
condition of zero macroscopic displacement field : 



no 



dV, 



/3 



(7) 



Introducing in Eq. ^ the relationship between field £ and polarization V , deduced from Eq. under the condition 
of vanishing displacement field, Born and Callen charges can be related to each other thanks to the knowledge of the 
optical dielectric tensor e°° : 
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SO that finally 
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For the case of isotropic materials, we recover the well known equality : Zk*""^' = e'^.Z^''^^ Even if they are both 
related to the change of polarization induced by an atomic displacement. Born and Callen charges appear as two 
distinct quantities and will be significantly different in materials where e°° is different from unity. 

Basically, an infinite number of charges could be defined corresponding to different boundary conditions, relating V 
and £. One of them is the Szigeti charge - Z*^''^^ -, defined as the change of polarization under the condition 

of vanishing local field, fjoc, at the atomic site: 



*(L) 



7*(S} 



no 



(10) 



£laa=0 



Contrary to Born and Callen charges, Z*'"^' was sometimes considered as a mo del- dependent concept in the sense 
that the local field is not observable as the macroscopic field but require some assumptions to be estimated. In the 
particular case of an isotropic material, the condition of vanishing local field can be written as follows: 



An 

£ioc=£+ -^V = 



Introducing this condition in Eq. (||) : 
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SO that we find : 
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In calculations of the dynamical properties of crystals, the contribution from the long-range Coulombic interaction 
to the atomic forces is usually restricted to dipolar forces and is included through a term : — 8. From 

Eq. ([T^), it can be checked that this force can be alternatively written in terms of local quantities : — Z*J-^^ Niac- 
in shell-model calculations, this second formulation is usually preferred. Indeed, from its definition, Z*^^^ only includes 
the effects of charge redistribution resulting from short-range interactions and it is therefore conveniently assimilated 
to the static charge [^5 22 1. 

From the previous discussion, it appears that the amplitude of the dynamical charge in macroscopic bodies is 
sensitive to the condition imposed on the macroscopic electric field. Considering finite clusters of increasing size, 
we deduce that the amplitude of the dynamical charge, reducing to the APT for a microscopic body, will tend to 
a different value when the macroscopic limit is taken, depending from the shape of the cluster. We investigate now 
this observation in more details, and provide a unified treatment of dynamical charges in periodic solids and clusters, 
sufficiently large for the macroscopic quantities (£, 7^, e°°, ...) to be defined. 

Following the well-known practice for the study of dielectric bodies , we consider that the cluster has a macro- 
scopic ellipsoidal shape. In this case, the macroscopic field within the cluster present the practical advantage to be 
homogeneous. In absence of any applied external field, it reduces to the depolarizing field related to the macroscopic 
polarization thanks to the depolarization coefficients Ua p8| . If we assume in what follows that the principal axes of 
the ellipsoid are aligned with the axes of coordinates, we have the following relationship : 



(14) 



where the geometry imposes : rii = \. Following the same procedure as previously, the dynamical charge Z*'^^ of 
a given atom k in an ellipsoid of volume U, can be written as : 
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and we have the general relationship : 
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In this expression, the presence of the depolarization coefhcients emphasizes the influence of the shape of the cluster 
on the amplitude of Z*^-^\ The above-mentioned sum rule on the depolarization coefficients forbid to impose the 
condition of zero electric or displacement fields simultaneously in the three directions. However, we have the following 
three interesting cases. 

First, we consider an extremely oblate ellipsoidal (slab- like) cluster and take the macroscopic limit. Along the z 
direction perpendicular to the surface, 1, while, along the two other directions, rix = riy ^ 0. The dynamical 

charge for the ellipsoid is therefore related to the Born effective charge through the following expression: 



^ K,,OLX 




(18) 



For uniaxial systems with no off-diagonal terms in the dielectric tensor, we note that the cluster charge along the 
direction perpendicular to the slab becomes identified with the Callen charge, while that in the slab plane reduces to 
the Born effective charge. 

Differently, for an extremely prolate ellipsoidal (needle-like) cluster aligned along the z direction (for which 
and Ux = ny 1/2), we have the following relationship: 
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(19) 



Here also, the charge along the z direction will reduce to the Born charge in uniaxial systems. 

Finally, for a spherical cluster, the symmetry imposes ni = ?i2 = = 1/3, so that = —AnV /i. For the case of 
an isotropic material, we recover therefore the condition of vanishing local field and becomes equivalent to ■ 
Therefore, we obtain the interesting result that in isotropic compounds, the Szigeti charge appears as a well-defined 
quantity and is simply the dynamical charge observed in a spherical cluster. 

To summarize, the concept of dynamical charge in macroscopic systems is not uniquely defined : it depends on the 
relationship between £ and V. In each case, the charge was however expressed in terms of two basic concepts, Z*^-^^ 
and e°°. In this Section, we focused on the term that includes the dielectric constant, and that describes the part 
of the electronic charge redistribution induced by the presence of a macroscopic field. In the next Section, we will 
discuss the physical processes responsible of the amplitude of Z*^'^\ 



B. Dynamical changes of orbital hybridizations 

During the seventies, a large variety of semi-empirical models were proposed to investigate the underlying physical 
processes driving the amplitude of dynamical charges. Without being exhaustive, let us mention the interesting 
treatments of Lucovsky, Martin and Burnstein |Q who decomposed Z**^^) in a local and a non-local contribution, of 
Lucovsky and White [ pO| discussing Z*^^) in connection with resonant bonding properties, or the bond charge model 
of Hiibner ||5l|. The most popular and sophisticated of these approaches remains however that of Harrison [|6 52-54 



within his bond orbital model (BOM). A similar theory was developed independently by Lannoo and Decarpigny [p5| . 

The BOM basically consists in a simplified tight-binding model, where the Hamiltonian is limited to the on-site and 
nearest-neighbour terms. The on-site elements are identified to free atom terms value, while the interatomic elements 
are taken as universal constants times a particular distance dependence. Among other things, these parameters 
determine the transfer of charge between the interacting atoms. As noted by Dick and Overhauser ||5^, the charge 
redistribution produced by the sensitivity of the overlap integrals on the atomic positions is at the origin of an 
"exchange charge polarization" . Similarly, in the Harrison model, the dependence of the parameter on the bond 
length are at the origin of dynamical transfer of charges and monitors the amplitude of Z*^"^^ that can become 
anomalously large as it is illustrated in the following examples. 

Let us first consider a diatomic molecule XY, composed of two open shell atoms, where Y has the largest elec- 
tronegativity. The interatomic distance is u and the dipole moment p{u). These observables allow us to identify a 
convenient static charge Z{u) = while the dynamical charge is defined as : 



Z*{u) 



dp{u) 
du 

= i?M + .2|M (20, 

In the last expression, Z* appears composed of two terms. The first one is simply the static charge. The second 
corresponds to an additional dynamical contribution: it originates in the transfer of charge produced by the modifica- 
tion of the interatomic distance. Within the BOM, this last contribution is associated to off-site orbital hybridization 
changes and is deduced from the universal dependence of the orbital interaction parameters on the bond length, as 
illustrated on a practical example in Appendix We deduce that the difference between Z{u) and Z*{u) will be 
large if Z(u) changes rapidly with u. It can even be non-negligible when dp{u)/du is small, if the charge is transferred 
on a large distance u. 

Moreover, this simple model naturally predicts anomalous amplitude of the dynamical charges, i.e. a value of Z*{u) 
not only larger than the static charge Z{u) but even larger than the "nominal" ionic charge. As the distance between 
X and Y is modified from to some u, the distance corresponding to a complete transfer of electrons from X to Y, 
the dipole moment evolves continuously from p(0) = (since there is no dipole for that case) to piu). Interestingly, 

Z* {u) du ^ [p{u) - p(0)] =uZ{u) (21) 
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so that: 



1 



Z*{u)du^ Z{u) (22) 



From the last relationship the mean value of Z*{u) from to It is equal to Z{u) (the "nominal" ionic charge). 
Consequently, if Z{u) changes with u, Z*{u) must be larger than Z{u) for some u between [0,?l]. The difference 
between Z*(u) and the nominal charge Z{u) is usually referred to as the anomalous contribution p7[ |. 

Switching now from a molecule to a linear chain ...-Y-X-Y-..., and displacing coherently the X atoms by du, shortened 
and elongated bonds will alternate all along the chain. For Harrison Q , the interaction parameters will be modified 
such that "the covalent energy increases in the shorted bond, making it less polar by transferring electron to the positive 
atom" . Inversely, electronic charge will be transferred to the negative atom in the elongated bond. These transfers 
of charge will propagate all along the chain, so that even if the net charge on the atom is not modified, a current of 
electrons will be associated to the atomic displacement. The direction of this electronic current is opposite to that of 
the displacement of positive atoms, so that it reinforces the change of polarization associated to this displacement and 
may generate an anomalously large dynamical charge. In our example, we have implicitely considered a truly periodic 
system under the condition of zero macroscopic electric field so that the associated dynamical charge is Z*^'^\ Under 
other conditions, the amplitude of the transfers of charge would be additionally influenced by the presence of the field 
as discussed in the previous Section. We note that, contrary to what was observed for the static charge, consequences 
of the covalence effects are to increase the amplitude of Z*^'^\ 

The previous model can finally be extended to three dimensional solids. For this case, however, the calculation 
of the dynamical contribution may become questionable when the identification of the charge transfers is restricted 
to some specific bonds [ p8| . As it will be discussed in Sections ^ and VI the Harrison model remains however a 
meaningful picture of practical interest to interpret more accurate results. 

Up to now, we focused on a "delocalized" model within which the electronic charge redistribution induced by an 
atomic displacement is visualized by transfer of charge induced by off -site changes of hybridization. In the past, 
various shell-models have however also been developed to investigate the dynamical properties of crystals. In these 
calculations, an accurate description of Z*^'^'> was mandatory in order to reproduce correctly the splitting between 
longitudinal and transverse optic modes in the vicinity of the F point. Contrary to the BOM, the shell-model is 
"local" and treats the charges within the Clausius-Mosotti limit. The previous dicussion in terms of a static and 
dynamical contribution to Z*'^'^^ remains valid. However, the dynamical contribution results there simply from the 
relative displacement of the shell charge as a whole with respect to the atom. It is attributed to the polarizability of 
the electrons in the local field at the atomic site. In the language of the BOM, such a displacement of the electronic 
cloud can be understood in terms of on-site changes of hybridizations. This approach contrasts with the model 
developped by Harrison but can also yield plausible Born effective charge amplitudes /citeGhosez95b. 

It must be emphasized that it is not possible to discriminate a priori between localized and delocalized models. 
Within the recent theory of polarization, it has been clarified that for the purpose of understanding polarization 
problems, "the true quantum mechanical electronic system can be considered as an effective classical system of 
quantized point charges, located at the centers of gravity associated with the occupied Wannier functions in each unit 



cell" |59| . Consequently, the correct description of the Born effective charges does not require to reproduce correctly 
all the features of the valence charge distribution but only the displacement of its Wannier center (see Ref. po|). 
As schematized in Fig. ^, antagonist models can reproduce a similar displacement of the Wannier center. In real 
materials, both local polarizability and transfers of charge do probably contribute to the charge reorganisation. It 
will be emphasized later how first-principles investigations can help to identify the dominant mechanism. 

In conclusion, this Section has shown that Z* is related to the static charge (see Eq. but does not restrict to 
it: Z* may also include an additional, important, dynamical contribution. Whatever the mechanism of the charge 
redistribution (localized or delocalized), the amplitude of the dynamical contribution cannot be estimated from the 
inspection of the electronic density alone. So, we partly disagree with Meister and Schwarz Q who suggested that all 
the charges including the GAPT are driven by the same underlying parameter. In what follows, based on firs-principles 
calculations, we illustrate on different examples that Z*^'^'> may become anomalously large and independent of the 
amplitude of the static charge Z. Moreover, two atoms with similar Z can also exhibit strongly different Z*^'^K 

IV. A FIRST-PRINCIPLES APPROACH 

A brief review of the most commonly used first-principles approaches for computing the Born effective charges has 
been reported recently [pOj . Going beyond semi-empirical approaches, ab initio techniques allow accurate prediction 
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of Z*^"^") in materials where its amplitude is not necessarily directly accessible from the experiment. Going further, the 
first-principles approaches are also offering a new opportunity to clarify the microscopic mechanism modulating the 
amplitude of Z*^'^^ without any preliminary hypothesis. As it will be illustrated in the following sections, it reveals 
particularly useful to understand the origin of anomalously large Z*^^^ in ABO3 compounds. 

The results of the present paper have been obtained in the framework of the density functional formalism p2| , |63[ . 
The exchange-correlation energy has been evaluated within the local density approximation, using a Pade paramctriza- 
tion of Ceperley- Alder homogeneous electron gas data |^ . Integrals over the Brillouin-zone were replaced by a 
sum on a mesh of 6 x 6 x 6 special k-points |6^,^^ (10 points in the irreducible Brillouin zone). The "all electron" 
potentials were replaced by the same ab initio, separable, extended norm-conserving pseudopotentials as in Ref. [ p2[ . 
The wavefunctions were expanded in plane waves up to a kinetic energy cutoff of 35 Hartree (about 4100 plane waves). 

The Born effective charges have been deduced from linear response calculations ||^, using a variational formu- 
lation |p|, |68|j69| to the density functional perturbation theory. The decomposition of individual contributions from 
separate groups of occupied bands has been performed following the scheme described in Ref. ||6C| ]. The parameters 
used for the calculations guarantee a convergency better than 0.5% on Z*''^'^^ as well as on each of its band-by-band 
contributions. 



V. THE CUBIC PHASE OF ABO3 COMPOUNDS 



A. General results 



Recently, the Born effective charge tensors of perovskite ABO3 compounds have been at the center of numerous 
investigations [pl]- |l5| , ^[|7CH 73|| . In the cubic phase, they are fully characterized by a set of four independent numbers. 
The charge tensor of the A and B atoms is isotropic owing to the local spherical symmetry at the atomic site. For 
oxygen, the local environment is tetragonal and two independent elements 0|| and 0_l must be considered, referring 
respectively to the change of polarization induced by an atomic displacement parallel and perpendicular to the B-0 
bond. In Table ||, we summarize the results obtained within different approaches for the cubic phase of BaTiOs. 

The first reliable estimation of Z*^"^^ in BaTiOa is probably due to Axe |10|, from empirical fitting to experimental 
mode oscillator strengths [ fz^ . In ABO3 compounds, Z*^'^^ cannot be determined unambiguously from the experiment. 
However, within some realistic hypothesis. Axe identified the independent elements of the effective charges of BaTiOs 
and already pointed out their two essential features. First, the oxygen charge tensor is highly anisotropic. Second, the 
charges on Ti and 0|| contain a large anomalous contribution (i.e. an additional charge with respect to the nominal 
ionic value of +2 for Ba, +4 for Ti and -2 for O). 

Both these characteristics are confirmed by the first-principles calculations. Our ab initio results, computed from 
linear response, are also in excellent agreement with those of Zhong et al. ||l^, obtained from finite differences of 
polarization. The charge neutrality sum rule, reflecting the numerical accuracy of our calculation, is fulfilled to within 
0.02. We note that the values of Z*^"^-* are also qualitatively reproduced from a shell-model calculation |^ and 
accurately predicted within the SCAD model ff^ . 

The anomalous amplitude of the dynamical charge, reported in this Section, is not a specific feature of BaTiOa. 
Similar computations of Z**--^^ were performed on different perovskite ABO3 compounds and they all reproduce the 
same characteristics than in BaTiOs. A non exhaustive list of these results is reported in Table III. We observe that 
the choice of the A atom has a rather limited influence on Zg ' and Z^^^ ' , which appear closely related to the B 
atom. While the nominal ionic charge of Ti and Zr is -1-4 in these compounds, the Born effective charge is between 
-1-7.08 and -1-7.56 for Ti, and approximately equal to -)-6.03 for Zr. For Nb, the ionic charge is -1-5, while the Born 
effective charge is between -t-9.11 and +9.37. Extending the investigations to WO3 in the reference cubic phase (defect 
perovskite structure), the ionic charge on W is equal to +6, while the Born effective charge reaches the much larger 
value of +12.51. For the class of perovskite ABOaCompounds, it can be checked that Z*^'^^ evolves quasi linearly with 
the nominal charge of the B atom . 

For materials containing Pb, the previous considerations remain valid but there are additional anomalies concerning 
Z^ and Zq^ . This feature is due to the more covalent bonding of lead with oxygen that was illustrated in 



Ref. [p5 84 1 . In what follows, we will not be concerned with these lead compounds. 
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B. Origin of the anomalous contributions 



The approximate reciprocity between Oy and B anomalous contributions suggests that they should originate in 
a global transfer of charge between B and O atoms as described in Section III B| . In Ref . , Harrison had in fact 
already suggested the existence of giant Born effective charges in perovskite materials. Being unaware of the earlier 
results of Axe, he had however no experimental evidence to corroborate his semi-empirical calculations. 

In Appendix we report results obtained within the Harrison model (it follows the method described for KCl 
in Ref. H, p. 334.). For SrTiOa, from the universal tight-binding parameters of Harrison, we get a value of —8.18 



for , making plausible the giant anomalous effective charges only by focusing on the dynamical changes of 

hybridization between occupied O 2s-0 2p states and the unoccupied metal d states. In BaTiOa, the hybridization 
between O 2p and Ti 3d orbitals is a well known feature, confirmed by various sources (experiments LCAO 
calculations ||3l|-|3^ and DFT results [^,^). In this context, it was therefore realistic to focus on O 2p - B d 
hybridization changes to explain intuitively large anomalous contributions |T^. 

At the opposite, it may therefore appear surprizing that model calculations which does not explicitly include transfer 
of charges are able to predict correctly the amplitude of the Born effective charges. For instance, in Table I, we observe 
that the values of Z*^-^^ are qualitatively reproduced by a shell- model calculation |6^. A similar agreement between 
ab initio and shell model results was highlighted for KNbOa In both cases, the calculation was performed within 
the "polarizability model" introduced by Bilz et al. [fzotl , which includes an anisotropic and non- linear polarizability of 
the O atoms. In the same spirit, at the level of the SCAD model, the Born effective charges are accurately reproduced 
while there is no explicit transfer of electrons between the different atomic sites. As discussed in Section HI, antagonist 
models can be invoked to explain the origin of anomalous contributions as soon as they globally reproduce a similar 
displacement of the Wannier center of the valence charge distribution. What appear as a macroscopic current along 
the Ti-0 chain within the BOM shows itself as an unusual polarizability of the oxygen atoms within the shell model. 

It was not possible to discriminate unambiguously between localized and delocalized model until Posternak et al. IT^ ] 
proposed a convincing proof of the crucial role of "off-site" hybridizations. Based on first-principles calculations, they 
demonstrated for KNbOs that the anomalous contribution to the charge of Nb and Oj| disappears if the hybridization 
between O 2p and Nb 4d orbitals is artificially suppressed. In a similar spirit, the inspection of the Wannier functions 
of BaTiOs and the analysis of their deformation under an atomic displacement reported by Marzari and Vanderbilt ]8l| ] 
confirm the predominant role played by the Ti 3d orbitals and the explanation introduced by Harrison. 

In the next Section, we propose a band- by-band decomposition of the Born effective charges . This technique 

appears as a tool of paramount importance to clarify the microscopic origin of anomalous contributions. Identifying 
the dynamical transfer of charges without any preliminary hypothesis on the orbitals that interact, it will allow to 
generalize the basic mechanism that was proposed by Harrison. 



VI. IDENTIFICATION OF DYNAMICAL CHANGES OF HYBRIDIZATION 

In ABO3 compounds, the electronic band structure is composed of well separated sets of bands. The hybridizations 
between the orbitals of the different atoms are relatively small and each band can be identified by the name of the 
main atomic orbital which contributes to this energy level in the solid. The Born effective charge is defined by the 
change of polarization associated to a specific atomic displacement. Our purpose will be to identify the contribution 
of each well separated set of bands to this change of polarization (ljj6^ . 



A. Reference configuration 

In Ref. we have described how band- by-band contributions to Z*^-^^ can be separated from each others. 

Moreover, it has been demonstrated that the contribution to Z^^^ from a single occupied band n can be interpreted 
as a change of polarization VLo^Vp = —2.iloAdp where Adp is the displacement in direction /? of the Wannier center 
of band n, induced by the unitary displacement of the sublattice of atoms k in direction a. 

In order to understand the origin of the displacement of the Wannier center of each band, it is helpful to define 
a reference configuration that corresponds to what we would expect in a purely ionic material. In such fictitious 
material, each band would be composed of a single non-hybridized orbital and the Wannier center of each band would 
be centered on a given atom. In absence of any hybridization, when displacing a given sublattice of atoms k, the 
Wannier center of the bands centered on the moving atoms would remain centered on it, while the position of the 
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center of gravity of the other bands would remain unaffected. The contributions of these two kind of bands to 
would therefore be —2 and electrons, respectively. 

In the real material, the anomalous contribution of a particular band m to a given atom k is defined as the additional 
part with respect to the reference value expected in absence of any hybridization: it reflects how the center of the 
Wannier center of band m is displaced relatively to the atoms when the sublattice n moves |Q . Considering each band 
as a combinaison of atomic orbitals, such a displacement of the Wannier center of a band with respect to its reference 
position must be attributed to hybridization effects: it is associated to the admixture of a new orbital character to 
the band. When the orbitals which interact are located on different atoms ("off-site" hybridization), the dynamical 
changes of hybridization can be visualized as transfers of charge. If the interacting orbitals are on the same atom 
("on-site" hybridization), the mechanism much looks like a polarizability. 

Rigorously, our band-by-band decomposition is performed within DFT and formally only concerns the Kohn-Sham 
particles. It seems however that the results are rather independent of the one-particle scheme used for the 
calculation so that the results presented here should give a good insight on the physics of the ABO3 compounds. 



B. BaTiOa 



Let us first applied the band-by-band decomposition to barium titanate. The band structure of BaTiOa is presented 
in Fig. ^. Results of the decomposition of Z*'''^'^ in the theoretical cubic structure of BaTiOa are reported in Table IV. 
The first line {Z^) brings together the charge of the nucleus and core electrons included in the pseudopotential. The 
other contributions come from the different valence electron levels. The sum of the band-by-band contributions on one 
atom is equal to its global effective charge while the sum of the contribution to a particular band from the different 
atoms is equal to —2 (within the accuracy of the calculation), the occupancy of this band. 

Focusing first on the titanium charge, we observe that the Ti 3s contribution (—2.03) is close to —2, confirming that 
the Ti 3s electrons follow the Ti atom when moving, independently from the change of its surrounding. This result 
a posteriori ]v&t\?ies the inclusion of deeper electronic levels as part of the ionic pseudopotentials. At the opposite, it 
is shown that the giant anomalous charge of titanium essentially comes from the O 2p bands (+2.86). It corresponds 
to a displacement of the Wannier center of the O 2p bands in opposite direction to the displacement of the Ti atom. 
This observation is in perfect agreement with the Harrison model: it can be understood by dynamical changes of 
hybridization between O 2p and Ti 3d orbitals, producing a transfer of electron from O to Ti when the Ti-0 distance 
shortens. This explanation was confirmed recently from the inspection of the O 2p Wannier functions Beyond 
the previous observations, we note however that there are also small anomalous charges from the Ti 3p, O 2s and Ba 
5p bands. These contributions are not negligible. The positive anomalous charges correspond to a displacement of 
the center of the Wannier function of the O and Ba bands in the direction of the closest Ti when this atom has moved. 



Some of these features go beyond the Harrison model, within which anomalous contributions to Zj^^ in Table IV 
would be restricted to the O 2p and O 2s bands. They suggest other kind of hybridization changes, that will be now 
more explicitly investigated. 

Focusing on barium, the global anomalous effective charge (-1-0.77) is small compared to that of Ti and this feature 
was first attributed to its more ionic character |13|. This ionicity is inherent to the Harrison model [|| and was 
confirmed in some ab initio studies Surprisingly, our decomposition reveals however that the anomalous 

charges of the O 2s (+0.73) and O 2p (+1.50) bands are not small at all. They are nevertheless roughly compensated 
by another Ba 5s (+0.11)and Ba 5p (+1.38) anomalous contributions. This result suggests that there are dynamical 
changes of hybridization between Ba and O orbitals as it was the case between O and Ti, except that the mechanism 
is here restricted to occupied states. This basically corresponds to a unitary transform within the suhspace of the 
occupied states which is unable to displace the global Wannier center of the valence charge. Our result so supports 
the hybridization of Ba orbitals, in agreement with experiment [^,^, LCAO calculations |3^,Q and DFT js^ j 
computations. Similar compensating contributions were recently observed in ZnO which has conventional Born 
effective charges |83) and in a series of alkaline-earth oxides . 

We note that a confusion sometimes appears that should be removed: the amplitude of the anomalous contributions 
to Z*^'^'^ is not related to the amplitude of the hybridizations but to the rate of change of these hybridizations under 
atomic displacements. It is clear that, in BaTiOa, the Ba 5p contribution to the O 2p bands is smaller than the 
contribution from the Ti 3d orbitals |^,^. However, the high sensitivity of this relatively weak covalent character 
under atomic positions is sufficient to produce large band by band anomalous contributions to Z*^^^. From that 
point of view, the Born effective charge appears therefore as a sensitive tool to identify the presence of even small 
hybridizations. 
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Finally, concerning the oxygen, even if 0|| and Oj_ are defined respectively for a displacement of O in the Ti and 

*(T) *(T) *fT) *fT) rH — li 

Ba direction, it seems only qualitative to associate Zq^^ with Zj^y and Zq^ with Zg^ as suggested in Ref. [ |13[ . 

The O 2p anomalous contributions to Ti and 0|| do not exactly compensate. Moreover, O 2p contribution to Z^^^ 
does not come from 0_l only but has equivalent contributions from 0||. This seems to confirm the idea of Bennetto 
and Vanderbilt that in 3D materials, transfers of charges are not necessarily restricted to a particular bond, but 
is a rather complex mechanism that must be treated as a whole. 

To summarize, our study has clarified the mixed ionic-covalent character of BaTiOs: it clearly establishes that 
the covalent character is not restricted to the Ti-0 bond but also partly concerns the Ba atom. Moreover, it leads 
to a more general issue. It illustrates that the presence of a large anomalous charge requires a modification of the 
interactions between occupied and unoccupied electronic state. The contributions originating from the change of the 
interactions between two occupied states correspond to unitary transforms within the subspace of the valence charge : 
they compensate, and do not modify the global value of Z*^'^\ 



C. SrTiOs 

The same analysis is now performed on SrTiOs. Its band structure (Fig. ^ is very similar to that of BaTiOa, 
except that the Ti 3p and Sr 4s bands are energetically very close to each others. Consequently, they strongly mix 
and it should be relatively meaningless to separate their respective contributions. The Sr 4p and O 2s states are in 
the same energy region but can be separated, contrary to what was observed in a study of SrO Q . 

The result of the decomposition is very similar (Table M) to that reported for BaTiOa. There is still a giant 
contribution to Zrp^ from the O 2p bands. On the other hand, while the Ba 5p bands were approximately centered 
between O 2s and O 2p bands in BaTiOs, the Sr 4p electrons are closer to the O 2s bands and mainly hybridize with 
them in SrTiOa. This phenomenon produces large but compensating contributions from Sr 4p and O 2s bands to 
Zg\. . Such an evolution is in agreement with the picture that anomalous contributions originate from off-site orbital 
hybridization changes. 



D. Other examples 

From the two previous results that concern two very similar materials, it might be suggested that not only the 
dynamical hybridization of the valence bands with unoccupied d-statcs but also the particular cubic perovskite struc- 
ture of ABO3 compounds plays a major role in determining Z*^-^^. For instance, the anomalous charge could partly 
originate in the local fields at the atomic sites, known to be anomalously large in this structure JSTt 

It is interesting to observe that anomalous charges are not restricted to perovskite solids but were also detected 
in a series of alkaline-earth oxides of rocksalt structure (CaO, SrO, BaO) jS^JS^ or even AI2RU all examples 

where the unoccupied d-states seem to play a major role. Interestingly, two materials belonging to the same structure 
can present completely different charges. This was illustrated for the case of Ti02 rutile and Si02 stishovite pl| , p2[ : 
while relatively conventional charges were observed on Si (-1-4.15) and O (-2.46) along the Si-0 bond in stishovite, 
giant effective charges, similar to those of BaTiOa, were obtained on Ti (-1-7.33) and O (-4.98) along the Ti-0 bond in 
rutile. Similarly, no anomalous charge was reported for MgO {Zq'^'^ = —2.07), presenting the same rocksalt structure 
than BaO {Zq'^'' = —2.80) In the same spirit, the same atom in different environments can present similar 

dynamical charge, as illustrated for Z*f ^ in BaTiOa and Ti02 H, or for Z*^^^ in BaZrOa (l|] and Zr02 H]. Also, 
in the family of ABO3 compounds, giant effective charges are observed on Ti in CaTiOa (Z*f' = 7.08, [||) but not 
on Si in CaSiOg (Z*f ^ = 4.00, ||). 

We observe that the presence of partly hybridized d-states seems the only common feature between the materials 
presenting giant anomalous effective charges, listed up to date. This feature finds a basic justification within the BOM 
of Harrison: the interaction parameters involving d-states are indeed much more sensitive to the interatomic distance 
than those involving, for example, s and p orbitals They will therefore be associated to larger dynamical transfer 
of charge and will generate higher Z*'-'^\ 
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VII. SENSITIVITY OF Z*'^' TO STRUCTURAL FEATURES 



In the litterature, calculations of Z*^'^^ essentially focused on the cubic phase of ABO3 compounds 1 11 15 6l] , |70| , |7l|| . 
On the basis of an early study of KNbOa ||ll[], it was concluded that the Born effective charges are independent of 
the ionic ferroelectric displacements (i.e. they remain similar in the different phases). Another investigation in the 
tetragonal phase of KNbOs and PbTiOa |l^, seemed to confirm that Z**^"^' are quite insensitive to structural details. 

These results were surprising if we remember that anomalous contributions to Z*'-^) are closely related to orbital 
hybridizations, these in turn, well known to be strongly affected by the phase transitions |^,|^. We will see in this 
Section that, contrary to what was first expected. Born effective charges in BaTiOa are strongly dependent of the 
structural features. 

We first investigate the sensitivity of the Born effective charges to the ferroelectric atomic displacements . For 
that purpose, we compute Z*^^-' in the three ferroelectric phases at the experimental unit cell parameters, with relaxed 
atomic positions as reported in Ref . . Table VI summarizes the results for a cartesian set of axis where the z-axis 
points in the ferroelectric direction. The Ba and Ti charge tensors are diagonal in each phase for this particular 
choice. In the case of O, we note the presence of a small asymmetric contribution for the lowest symmetry phases. 
The eigenvalues of the symmetric part of the tensor are also reported. In each phase, the eigenvector associated to 
the highest eigenvalue of O approximately points in the Ti-0 direction and allows to identify the highest contribution 
as 0||. The other eigenvalues can be referred to as Oj_, by analogy with the cubic phase. 

Although the charges of Ba and 0_l remain globally unchanged in the 4 phases, strong modifications are observed 
for Ti and 0||: for example, changing the Ti position by 0.076A (2% of the unit cell length) when going from the 

cubic to the rhombohedral phase, reduces the anomalous part of Z^^^^ by more than 50% along the ferroelectric axis 
(Table |Vl|). Equivalent evolutions are observed in the other ferroelectric phases. Similar changes were detected in 
KNb03~||l. 

In the isotropic cubic structure, Harrison had explained the large value of Z*^^-' in terms of the Ti-0 bond length. 
For the anisotropic ferroelectric phases, it should be intuitively expected that the shortest Ti-0 distance dmm in 
the structure will dominate the bonding properties. It is therefore tempting to transpose the Harrison model to 
understand the evolution of Z*^^^ in terms of the distance dmin- In Fig- H, the amplitude of Z^^^^ in the direction 
of the shortest Ti-0 bond length of each phase is plotted with respect to dmin- I similar graph can be obtained for 
O. For the different phases, at the experimental lattice parameters, we observe that the anomalous parts evolve quasi 
linearly with dmin- 

Independently from the previous calculations, we also investigated the evolution of Z*^"^^ under isotropic pressure 
(Table VH). In contrast with the changes observed with respect to the atomic displacements, the charge appears 
essentially insensitive to isotropic compression. In particular, in the compressed cubic cell at 3.67 A where the Ti-0 
distance is the same that the shortest Ti-0 bond length in the tetragonal structure , Z^^^"^ remains very close to 
its value at the optimized volume. This new element clearly invalidates the expected dependence from Z*^^-* to dmin- 

The fundamental difference between the cubic and tetragonal structures lies in the fact that in the cubic phase 
every Ti-0 distance is equal to the others, while in the tetragonal phase, along the ferroelectric axis, a short Ti-0 
bond length (dmin) is followed by a larger one (dmax) which breaks the Ti-0 chain in this direction. In order to verify 
that it is not this large Ti-0 distance which, alternatively to dmin, is sufficient to inhibit the giant current associated 
to the anomalous charges, we also performed a calculation in an expanded cubic phase where Oo = 2.dmax- we observe 
however that the Ti charge is even larger than in the optimized cubic phase. 

We conclude from the previous investigations that the amplitude of Z*^'^^ in BaTiOz is not dependent on a particular 
interatomic distance (dmin, dmax) but is more critically affected by the anisotropy of the Ti environment along the Ti- 
O chains. In agreement with this picture, Wang et al. j?^ reported recently an insensitivity of Z**^-^) to a tetragonal 
macroscopic strain in KNbOs. Also, the charges reported by Bellaiche et al. in mixed a compound as PZT, where 
the ionic environment becomes anisotropic, seem to confirm our results. 

A band by band decomposition of Z^^^^ (Table VIII ) points out that the difference between the cubic and tetragonal 
phases is essentially localized at the level of the O 2p bands (-1-1.48 instead of -1-2.86) while the other contributions 
remain very similar. This suggests an intuitive explanation. In the cubic phase the O 2p electrons are widely delocal- 
ized and dynamical transfer of charge can propagate along the Ti-0 chain as suggested by Harrison. In the tetragonal 
phase, the Ti-0 chain behaves as a sequence of Ti-0 dimers for which the electrons are less polarizable. This smaller 
polarizability is confirmed by a similar reduction of the optical dielectric constant along the ferroelectric direction. 
This analysis seems plausible from the Wannier function analysis reported recently by Marzari and Vanderbilt 

Finally, let us mention that if the evolution of Z*^-^) is relatively weak under isotropic pressure, it would be wrong to 
consider that the dynamical properties of BaTiOs are insensitive to the volume: small changes are observed that are 
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of the same order of magnitude than for other compounds hke SiC |£8|,^ . The direction of the evolution is however 
different. Moreover, the evolution of the different charges is even not identical: while the absolute value of Zg^ and 
Z'qJ^ decreases with increasing volume, the inverse behaviour is observed for Z'^J^ and Z^^\ 

Here also, the band by band decomposition (Table fx\) reveals some hidden features. In the compressed cubic 
phase, the anomalous part of the Ba 5p , Ba 5s and Ti 3p bands are 50% larger than in the optimized cubic cell. 
This suggests an evolution of the interactions between occupied orbitals that is coherent with the modification of the 
interatomic short-range forces observed independently At the opposite, in our expanded cubic phase, most of the 
anomalous contributions to Zg^ and Zj,\ have disappeared in agreement with the picture of a more ionic material. 
The O 2p contribution, is the only one that remains surprisingly large. Comparing to the value obtained for the cubic 
phase at the experimental volume, its evolution was even more important than the linear dependence upon the bond 
length, expected from the Harrison model. 



VIII. SPONTANEOUS POLARIZATION 



The spontaneous polarization (P^) of the ferroelectric phases can be determined by integrating the change of polar- 
ization along the path of atomic displacement from the paraelectric cubic phase (taken as reference) to the considered 
ferroelectric structure. If the effective charges were roughly constant, this integration should be approximated by: 

p^-o. = 7^- E ^^-'P (23) 

However, we have seen, in the previous Section, that the Born effective charges are strongly affected by the atomic 
displacements. It is therefore important to investigate their evolution all along the path of atomic displacements from 
one structure to the other. 

We performed the calculation for a transformation from the cubic to the rhombohedral structure. The rhombohedral 
macroscopic strain is very small and was neglected |lO0|| : our calculation was performed by displacing the atoms to 
their theoretically optimized position in rhombohedral symmetry, when keeping the cubic lattice parameters. The 
result is reported in Figure |6|, for Zj}- ' along the ferroelectric direction. A similar curve can be obtained for ' . 

We observe that the evolution of Z*^^^ is approximately quadratic close to the cubic phase. However, it becomes 
rapidly linear, and remains linear for displacements even larger than those associated to the ferroelectric distortion. 

Expecting a similar evolution of the dynamical charges for the tetragonal and orthorhombic displacements, an 
estimation of the spontaneous polarization in the ferroelectric phases can be found when using Eq. ( p3|) with a mean 
effective charge determined from its value in both phases. Using a mean charge estimated from the values in the para- 
and ferro-electric phases, we obtain the spontaneous polarization s pre sented in Table 



Our results are only in relative agreement with the experiment |101,r^ and suggest different comments. Firstly, we 



would like to mention that part of the discrepancy must be assigned to the theoretical overestimation of the computed 



ferroelectric displacements, discussed in Ref. [g5| : when using the experimental displacements of Ref. 1102], we 
recover a better estimation of Pg as in Ref. ||l^ . The dispersion of X-rays diffraction data makes however difficult the 
exact identification of the ferroelectric displacements. Secondly, another part of the error could be due to the lack of 



polarization dependence of the LDA |102]. Finally, we note that there is also some uncertainty on the experimental 
value of Pt. 



IX. CONCLUSIONS 



In this paper, we first analyzed the links between different definitions of atomic charge. We have shown that, 
contrary to the static definitions, dynamical effective charges also depend on the electronic charge reorganisation 
induced by an atomic displacement. The amplitude of this dynamical contribution is monitored not only by the 
bonding with the other atoms but also, for large systems, by the condition imposed on the macroscopic electric field. 
A unified treatment of the concept of dynamical charge in molecules, large clusters, and truly periodic systems has 
been presented, in which the Born effective charge and the optical dielectric constant appear as the two fundamental 
quantities. The microscopic origin of the dynamical contribution has been discussed in terms of local polarizability 
and delocalized transfers of electrons. 
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Based on various first-principles results, we have then emphasized that the Born effective charges are anomalously 
large in the family of ABO3 compounds: their amplitude can reach more than twice that of the nominal ionic charges. 
This feature was explained in terms of interatomic transfers of charge, produced by "off-site" dynamical changes 
of hybridization. For BaTiOs and SrTiOa, we have brought to light complex dynamical changes of hybridization, 
concerning not only Ti and O but also Ba and Sr orbitals. The hybridizations restricted to occupied states generate 
however compensating anomalous contributions so that, finally, the total value of Z*'^-^-' is essentially affected by 
dynamical changes of hybridization between O 2p and Ti 3d orbitals. 

As a more general issue, it appears that the existence of partial hybridizations between occupied and unoccupied 
states is an important feature for candidate to large anomalous Born effective charges. Moreover, the dynamical 
transfers of charge are expected to be larger when such a hybridization involve d states, for which the interactions 
parameters with other orbitals are particularly sensitive to the interatomic distance. 

Investigating the evolution of Z*'^'^^ to the structural features, we have shown that they are strongly affected by 
the ferroelectric atomic displacements and much less sensitive to isotropic pressure. The results have clarified that 
the amplitude of Z*^"^^ is not monitored by a particular interatomic distance but is dependent on the anisotropy of 
the Ti environment along the Ti-0 chains. 

Finally, the effective charges were used to estimate the spontaneous polarization in the ferroelectric phases of 
BaTiOa. For that purpose, their evolution was investigated all along the path of atomic displacement from the cubic 
to the rhombohedral structure. 

All along this work, we only focused on the microscopic mechanisms that govern the amplitude of the Born effective 
charges. In independent studies, it was also emphasized that the anomalously large Born effective charges produce 
a giant LO-TO splitting in ABO3 compounds, specially for the ferroelectric phonon mode jljjl^. Moreover, it 
was demonstrated that this feature is associated to the existence of an anomalously large destabilizing dipole-dipole 
interaction, sufficient to compensate the stabilizing short-range forces and induce the ferroelectric instability . In 
materials where polar modes play a major role, the Born effective charge appears therefore also as a "key concept" 
to relate the electronic and structural properties |104|. 
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APPENDIX A: THE HARRISON MODEL 



In this Appendix, we briefly describe the bond-orbital model proposed by Harrison for the case of ABO3 com- 
pounds 1^]. In particular, we pay a particular attention to the definition and the calculation of static and dynamical 
charges within this model. Values are reported for SrTiOa. 

The bond orbital model of Harrison consists in a simplified tight-binding model, where the Hamiltonian is limited 
to the on-site and nearest-neighbour terms. Moreover, for ABO3 compounds, it is assumed that the A atom has no 
other function than to provide electrons to the system, and is fully ionized. The only considered interactions involve 
B and O atom orbitals. 

The model includes O 2s, O 2p and B d orbitals, interacting through Vsd,a, Vpd,a and Vpd.-K parameters. As the 
matrix elements Vsd,a and Vpd.a are nearly identical, it is suggested to construct two sp hybrids on the oxygen, 
\h± >= (|s > ±.\p >)/\/2. Each of these hybrid will have a large matrix element (V^^ ^ = \Vsd,a + Vpd,a\/ V^) coupfing 

it to the d state on one side and a negligible matrix element {Vi~^ ^ — \Vsd.ij — Vpd,a\/ V^) coupling it to the d state on 
the other. 
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1. Effective static charge 



In absence of orbital interactions, the static charges of SrTiOa would be of +2 on Sr, +4 on Ti and —2 on O. 
However, due to the Ti-0 orbital interactions, the transfer of electrons from Ti to O is not complete. Within the 
BOM, the effective static charge on O can be estimated as: 

Zo = -2 + T, + (Al) 

where 

("4" because there are 2 hybrids composed of 2 electrons that each interacts mainly with one Ti neighbour) and 

T^ = 8(j^^] (A3) 

("8" because there are 4 electrons that are each partly delocalizcd on the 2 Ti neighbours). 

For SrTiOa, from the parameters of Matheiss (Ref. p. 445), = 0.35 and T^r = 0.68 so that the effective static 
charge on the oxygen atom is Zq = 0.96. 



2. Born effective charge 

Following the discussion of Section III-B, the Born effective charge can be obtained by adding to the static charge 
the dynamical contribution originating in the dependence of the matrix elements V upon the interatomic bond length 
d. From Harrison's Solid State Table (see also Eq. (19-11) in Ref. |^), the previous matrix elements have all the same 
interatomic dependence: 

so that 

{V + SVf ^V^ + 2V SV + 0{SV'^) (A6) 
= V'^ + 2{^)V^^+0{5V^) (A7) 

When displacing the O atom along the Ti-0 direction, there will be an additional transfer of electron from O to 
the nearest Ti that is equal to 

5T.^2[2(^)(^^)^]^ (A8) 

<5T. = 4[2(^)(^^)2]^ (A9) 

On the other hand, there will be the same transfer of charge from the other neighbour Ti atom to O, so that the 
previous electrons are globally transferred on a distance 2d. The change of polarization associated to this transfer of 
charge is: 

dP, = -28( ^''^'^ f5d (AlO) 
Ed - e/i 

5P^ = -56(-5^i^)2<5c/ (All) 

Ed — Ep 

The associated dynamic contribution to the Born effective charge is (SP/ Sd) : 



15 



Sz<T) ^ „28(-5^)2 (A12) 

Cd — 

SZ<^^ = -56(-^^)2 (A13) 

The effective charge on the O atom for a displacement along the Ti-0 direction is therefore: 

Z;(p=Zo+5Z:(^)+<5Z:(^) (A14) 

For SrTiOa, from the parameters of Matheiss (Ref. |], p. 445), 8Z*}^^ = -2.45 and iZ*}'^^ = -4.76 so that 
Z*f ^ = -8.18. 
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FIG. 1. Ti-0 plane cut in the partial electronic density coming from the O 2p bands in the cubic phase of barium titanate. 




FIG. 2. Schematic representation of the two basic mechanisms that can explain the displacement of the Wannier center of a 
band under atomic displacement : (a) local polarizability, (b) interatomic transfers of charge. 
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FIG. 3. Kohn-Sham electronic band structure of BaTiOa. 




FIG. 4. Kohn-Sham electronic band structure of SrTiOs. 
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TABLE I. Static charges of BaTiOs in the cubic structure. 



ZBa 



Zq 



Reference 



Nominal 

Empirical models 



First-principles calculations 



+2 
+2.00 
+ 1.40 
+2.00 
+1.86 
+1.48 
+2.00 
+2.00 
+2.12 
+1.39 



+4 
+0.19 
+2.20 
+1.88 
+3.18 
+1.86 
+1.86 
+2.89 
+2.43 
+2.79 



-2 
-0.73 
-1.20 
-1.29 
-1.68 
-1.11 
-1.29 
-1.63 
-1.52 
-1.39 



Ref. ji 
Ref. [17 1 
Ref. [22 1 
Ref. [191 
Ref. [20 1 
Ref. [n| 
Ref. [23 1 
Ref. [M 
Ref. [~| 



TABLE II. Born effective charges of BaTiOs in the cubic structure. 



*(T) 



,{T) 



*(T) 



Reference 



Nominal 

Experiment 

Models 

First-principles 



(Shell model) 
(SCAD model ) 
(Linear response) 
(Berry phase) 



+2 
+2.9 
+1.63 
+2.9 
+2.77 
+2.75 



+4 
+6.7 
+7.51 
+7.3 
+7.25 
+7.16 



-2 
-2.4 
-2.71 
-2.2 
-2.15 
-2.11 



-2 
-4.8 
-3.72 
-5.8 
-5.71 
-5.69 



Ref. [10 1 
Ref. [n| 
Ref. [^1 
Present 
Ref. [|| 



TABLE III. Born effective charges of various ABO3 compounds in their cubic structure. 



AB03 




nominal 


2 


CaTiOa 


2.58 


SrTiOa 


2.56 




2.54 




2.55 




2.4 


BaZrOs 


2.73 


PbTiOs 


3.90 


PbZrOa 


3.92 


nominal 


1 


NaNbOg 


1.13 


KNbOs 


0.82 




1.14 




1.14 



*(T) 



*{T) 



Oil 



Reference 



4 
7.08 
7.26 
7.12 
7.56 
7.0 
6.03 
7.06 
5.85 



-2 
-5.65 
-5.73 
-5.66 
-5.92 
-5.8 
-4.74 
-5.83 
-4.81 



-2 
-2.00 
-2.15 
-2.00 
-2.12 
-1.8 
-2.01 
-2.56 
-2.48 



Ref. jl| 
Present 
Ref. jlS] 
Ref. |~] 
Ref. IK] 
Ref. i~] 
Ref. i~] 
Ref. |~1 



5 

9.11 
9.13 
9.23 
9.37 



-2 
-7.01 
-6.58 
-7.01 
-6.86 



-2 
-1.61 
-1.68 
-1.68 
-1.65 



Ref. lis) 
Ref. inj 
Ref. i~] 
Ref. l") 



nominal 
WO3 



6 

12.51 



-2 
-9.13 



-2 
-1.69 



Ref. 78 
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TABLE IV. Band by band decomposition of Z*^'^' in the optimized cubic phase of BaTiOs. The contributions have been 
separated into a reference value and an anomalous charge (see text). 











7(r) 


_L(JLcli 


z 

Ti 3s 

Ti 3p 
Ba 5s 
2s 
Ba 5p 
2p 


+10.00 
+ 0.01 
+ 0.02 

-2-0.11 
+ 0.73 

-6 - 1.38 
+ 1.50 


+12.00 
-2 - 0.03 

-6 - 0.22 
+ 0.05 
+ 0.23 
+ 0.36 
+ 2.86 


+6.00 
+ 0.00 

0-0.02 
+ 0.02 

-2 - 0.23 
+ 0.58 

-6 - 0.50 


+6.00 
+ 0.02 

+ 0.21 
+ 0.01 

-2 - 2.51 
0-0.13 

-6-3.31 


-2.00 

-6.03 
-2.01 
-6.01 
-5.99 
-17.95 


Ttjtal 


+ 2.77 


+7.25 


— 2.1 


-."').71 


-0.01 


TABLE V. 
separated into 


Band by band decomposition of Z**-^-* in the 
a reference value and an anomalous charge (i 


experimental cubic phase 
see text). 


of SrTiOs. The contributions have been 












Total 


Ti 3s 


+10.00 
+ 0.01 


+ 12.00 
-2-0.03 


+6.00 
+ 0.00 


+6.00 
+ 0.03 


+40 
-1.99 


Sr4s -1 
Ti3p J 
2s 

Sr 4p 
O 2p 


-2 + 0.02 


-6-0.18 


- 0.03 


+ 0.23 


-7.99 


+ 3.08 
-6-3.11 
+ 0.56 


+ 0.02 
+ 0.37 
+ 3.08 


-2-1.31 
+ 1.42 
-6-0.12 


-2 - 0.48 
0-0.10 
-6 - 3.41 


-6.00 
-6.00 
-18.01 


Total 


+2.56 


+7.26 


-2.15 


-5.73 


+0.01 
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TABLE VI. Born effective charges in the three ferroelectric phases of BaTiOs. Tensors are reported in cartesian coordinates, 
with the 2-axis along the ferroelectric direction. For Ba and Ti, the tensors are diagonal and only the principal elements 

are mentioned. For O, full tensors are reported. The eigenvalues of the symmetric part of ^*('^' arc mentioned in brackets; 
the eigenvector associated to the highest eigenvalue approximately points in the Ti direction. In the cubic phase, we had: 
Z'^P = 7.29, Z'Jp = 2.74, Z'^^ = -5.75 and Z^^ = -2.13. 



Tetragonal 



Orthorhombic 



Rhombohedral 



7*(T) 
-'Ba 

7*(T) 



7* (J 



(T) 



,(T) 



( +2.72 +2.72 +2.83 ) 

( +6.94 +6.94 +5.81 ) 

/ -1.99 \ 

j -1.99 I 

V -4.73 / 

[ -1.99 -1.99 -4.73 ] 

/ -2.14 \ 

I -5.53 

Vo -1.95/ 

[ -1.95 -2.14 -5.53 ] 

/ -5.53 \ 

j -2.14 

V -1.95 / 

[ -1.95 -2.14 -5.53 ] 



( +2.72 +2.81 +2.77 ) 

( +6.80 +6.43 +5.59 ) 

/ -2.04 \ 

I -3.63 +1.38 I 

V +1.57 -3.17 / 

[ -1.91 -2.04 -4.89 ] 

/ -2.04 \ 

-3.63 +1.38 

\ +1.57 -3.17 / 

[ -1.91 -2.04 -4.89 ] 

/ -5.44 \ 

-1.97 

V -2.01 / 

[ -1.97 -2.01 -5.44 ] 



(+2. 


.79 


+2. 


.79 


+2, 


(+6. 


.54 


+6 


.54 


+5, 




.54 


-0, 


.99 


+0, 


(1: 


.99 


-3. 


.68 


+1. 


V+o. 


.72 


+1. 


.25 


-2, 


f-i. 


.97 


-1. 


.98 


-5. 



-2, 


.54 


+0.99 


+0, 


+0, 


.99 


-3.68 


-1, 


+0. 


.72 


-1.25 


-2, 


-1, 


.97 


-1.98 


-5, 


-4, 


.25 





-1, 







-1.97 





-1.44 -2.78 / 
-1.97 -1.98 -5.05 ] 



TABLE VII. Evolution of the Born effective charges of BaTiOa under isotropic pressure in the cubic phase. 



tto (A) 




7*(T) 


7*(T) 


Z<^^ 


3.67 


+2.95 


+7.23 


-2.28 


-5.61 


3.94 


+2.77 


+7.25 


-2.15 


-5.71 


4.00 


+2.74 


+7.29 


-2.13 


-5.75 


4.40 


+2.60 


+7.78 


-2.03 


-6.31 
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TABLE VIII. Band by band decomposition of Zj,- in different structure of BaTiOs. 
into a reference value and an anomalous charge (see text). 



The contributions have been separated 



Band 


(cubic - 3.67 A) 


(cubic - 


- 3.94 A) 


(tetragonal - exp) 


7*{T) 

(cubic - 4.40 A) 


Ti 3s 
Ti 3p 
Ba 5s 
O 2s 
Ba 5p 
O 2p 


^ .uu 

-2-0.07 
-6 - 0.43 
+ 0.09 
+ 0.27 
+ 0.64 
+ 2.73 




-2 - 0.03 
-6-0.22 
+ 0.05 
+ 0.23 
+ 0.36 
+ 2.86 


4-1 9 on 
-2-0.05 
-6-0.26 
+ 0.05 
+ 0.25 
+ 0.34 
+ 1.48 


4-1 9 nn 
-2 + 0.01 
-6-0.07 
+ 0.02 
+ 0.19 
+ 0.13 
+ 3.50 


Total 


+7.23 




+7.25 


+5.81 


1 T 'TO 
+ 1 . 10 


TABLE IX. Band by band decomposition of Z^ J^ 
structure. The contributions have been separated into a 


in the optimized cubic phase of BaTiOs and in an 
reference value and an anomalous charge (see text). 


expanded cubic 


Band 




(cubic 


- 3.94 A) 




(cubic - 4.40 A) 


Z^ 

Ti 3s 
Ti 3p 
Ba 5s 
O 2s 
Ba 5p 
O 2p 






+10.00 
+ 0.01 
+ 0.01 

-2 - 0.11 
+ 0.73 

-6- 1.38 
+ 1.50 




+ 10.00 
0-0.01 
+ 0.01 

-2 + 0.00 
+ 0.37 

-6 - 0.44 
+ 0.66 


Total 






+2.77 




+2.59 



TABLE X. Spontaneous polarization in the three ferroelectric phases of BaTiOs in /iC/cm^. The results were de- 
duced from Eq. (G3) when using either Z*'^' from the cubic phase (Cubic) or a mean charge (Mean) defined as 



= 0.68 X Z: 



XT) 



+ 0.32 X Z 



ferroelectric displacements. 



ferroJ 



Results are reported for the experimental (Exp) and theoretical (Theo) atomic 



Positions 



Tetragonal 



Orthorhombic 



Rhombohedral 



Reference 



Cubic 
Cubic 
Mean 
Mean 



Exp 
Theo 
Theo 
Exp 



26.3 
30 
36.35 
34.02 
28.64 



30.7 
26 
42.78 
39.68 
36.11 



33.5 
44 
43.30 
40.17 



Exp. mm 

Ref. |3f 

Present 
Present 
Present 
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